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We present combined experimental and theoretical results on strong-field ionization of oriented 
carbonyl-sulphide molecules by circularly-polarized laser pulses. The obtained molecular frame 
photoelectron angular distributions show pronounced asymmetries perpendicular to the direction of 
the molecular electric dipole moment. These findings are explained by a tunneling model invoking 
the laser-induced Stark shifts associated with the dipoles and polarizabilities of the molecule and 
its unrelaxed cation. The focus of the present article is to understand the strong-field ionization 
of one-dimensionally-oriented polar molecules, in particular asymmetries in the emission direction 
of the photoelectrons. In the following article (Phys. Rev. A 83, 023406 (2011)) the focus is to 
understand strong-field ionization from three-dimensionally-oriented asymmetric top molecules, in 
particular the suppression of electron emission in nodal planes of molecular orbitals. 

PACS numbers: 33.80.Rv, 33.80.Eh, 42.50.Hz, 37.20. +j, 37.10.Vz 



I. INTRODUCTION 

When molecules are exposed to intense femtosecond 
laser pulses they ionize. If the ionizing laser pulse is lin- 
early polarized the electron can be steered back to rescat- 
ter on the ion left behind, thereby initiating phenom- 
ena such as high harmonic generation, above threshold 
ionization, and double ionization (see, e.g., Refs.fT^). 
Being the event that initiates these central strong-field 
processes, ionization has attracted special attention and 
much effort has gone into describing and understanding 
it. Because molecules are not spherically symmetric the 
ionization probability and the emission direction of the 
electron depends on the relative orientation between the 
molecule and the polarization vector of the laser pulse 
[5HT8] . Knowledge of this orientational dependence is im- 
portant for understanding, optimizing, or utilizing sub- 
sequent strong- field processes O |4] . 

The ability to align molecules, i.e., to confine one or 
more molecular axes along space-fixed axes, has over the 
past few years provided a valuable tool to experimentally 
explore the orientational dependence of strong-field ion- 
ization, and consequently, opened up for comparing the- 
oretical and experimental results [71 [TOl [11] |14l [15] |19] . 
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The majority of studies have focused on nonpolar lin- 
ear molecules where nonadiabatic alignment, by a lin- 
early polarized laser pulse, provides a convenient way to 
prepare a sample of 1-dimensionally aligned molecules 
[20l [21 . Most molecules are, however, polar, i.e., they 
do not exhibit inversion symmetry. The experimental 
investigation of such systems requires that not only the 
axes of the molecule are confined, the permanent electric 
dipole moment must also point in a particular direction. 
Thus, the molecule should be oriented in addition to be- 
ing aligned. Orientation can be achieved by static elec- 
tric field methods such as hexapole focusing [22[ |23] and 
brute- force orientation [24 , by optical methods based on 
two-color laser fields [25] [26] , or by combined laser and 
static electric field methods [27H3Q] . In the present work 
we employ the method relying on mixed laser and static 
electric fields since it provides very high degrees of align- 
ment and orientation. 

In detail we present a combined experimental and 
theoretical study of single ionization of a polar linear 
molecule, carbonyl-sulfide (OCS), by near infrared 30 
femtosecond laser pulses. First results were recently pre- 
sented elsewhere [31] . Unlike studies aimed at recollision 
phenomena, e.g., high-order-harmonic generation, circu- 
larly polarized pulses are employed. By doing so the 
strong-field dynamics is simplified since the circularly po- 
larized field drives the electrons away from the parent 
molecule and thus turns off recollision. Our studies fo- 
cus on the photoelectron angular distributions (PADs) 
from single ionization. When the OCS molecules are 



tightly aligned and oriented pronounced asymmetries are 
observed in the experimental PADs perpendicular to the 
fixed molecular axis. The asymmetries are absent for 
randomly oriented molecules. Our theoretical analysis, 
based on a modified tunneling theory, rationalizes the ex- 
perimental findings and shows that the observed asym- 
metries are determined by the difference in ionization 
probability between the cases when the circularly polar- 
ized field points in the same and in the opposite direction 
as the permanent dipole moment [32^ . Notably, the PADs 
refiect the permanent dipole moment and the polarizabil- 
ity of the active molecular orbital, which again, in the 
case of the highest-occupied molecular orbital (HOMO), 
may be related to the dipole moments and polarizabil- 
ities of the neutral molecule as well as of its unrelaxed 
cation. The calculated results are exponentially sensi- 
tive to (temporal) changes in these quantities and hence 
point to the extension to time-resolved measurements of 
valence electron dynamics using pump-probe settings. 

The paper is organized as follows: In Sec. II the ex- 
perimental technique is discussed and in Sec. Ill, the ex- 
perimental results are presented. Section IV presents the 
theory and in Sec. V the theoretical predictions are com- 
pared to the experiment. Conclusions are given in the 
last section. Appendix A provides a summary of the 
molecular properties of OCS and its cation relavant for 
the present study. 



II. EXPERIMENTAL METHOD 

The experimental setup is described in detail else- 
where ^33j so the discussion here is brief. A gas mixture 
of ~10 mbar carbonylsulfide (OCS) and 10 bar of Ne is 
expanded supersonically into vacuum through an Even- 
Lavie valve [3^ .35^ forming a pulsed molecular beam. 
The molecular beam is skimmed twice before entering a 
15-cm-long electrostatic defiector that spatially disperses 
the molecular beam in the vertical direction according to 
the quantum states populated [29l |30] . After exiting the 
defiector the molecular beam is crossed at 90° by two fo- 
cused laser beams, one to align and orient the molecules 
and one to induce ionization. The experiments described 
here are conducted on the most defiected molecules, i.e., 
a subset of molecules selected in the lowest lying rota- 
tional quantum states [29, 30 . 

The alignment beam originates from an injec- 
tion seeded Q-switched Nd:YAG laser (20 Hz, 
'?"FWHM = 10 ns, A = 1064 nm). The spotsize, {uj^^, 
in the focus at the crossing with the molecular beam is 
34 /im, yielding a peak intensity of ~ 8x10^^ W/cm^. 
The ionization laser beam, termed the probe beam, 
originates from a pulsed femtosecond Ti-Sapphire system 
(1 kHz, A = 800 nm) externally compressed to 30 fs 
(FWHM) and focused to cJq^^ ® = 21 /im resulting in a 
peak intensity of ~ 5.4x10^^ W/cm^. The probe pulse 
is electronically synchronized to the peak of the YAG 
pulse, where the degree of alignment is highest. 



The ions or electrons produced by the probe pulses are 
extracted with a weak static electric field in a velocity 
map imaging (VMI) geometry and projected onto a two 
dimensional detector consisting of a micro channel plate 
(MCP) detector backed by a phosphor screen. The ion 
or electron images on the phosphor screen are recorded 
by a CCD camera and the coordinates of each individual 
particle hit are determined. 



III. EXPERIMENTAL RESULTS 

A. Alignment and orientation 

The target of adiabatically oriented molecules is ob- 
tained by the combined action of the ac electric field from 
the YAG pulse and the weak static electric field present 
in the VMI spectrometer as shown in previous studies 
[29) |3Q] . The strongest orientation is obtained when the 
YAG pulse is polarized along static electric field, Fgtat 
and this is the geometry used in the photoelectron studies 



presented in Sec. IIIB To characterize orientation using 
2D ion imaging, as we do here, it is, however, necessary 
to rotate the molecules away from Fgtat as shown in Fig. 
Ilia). Our method to characterize alignment and orienta- 
tion is based on Coulomb exploding the molecules with 
an intense probe pulse and subsequently recording the ve- 
locities of the recoiling ions by the 2D imaging detector. 
This method does not work well for molecules aligned 
along Fgtat because all recoiling ions tend to collapse in 
the center of the detector. Therefore, the measurements 
of alignment and orientation are performed for molecules 
with their molecular axis rotated 45 degrees away from 

-f stat • 

Figure [l[b) and (c) display S+ ion images from 
Coulomb explosion of the OCS molecules with the probe 
pulse linearly polarized vertically and the YAG pulse po- 
larized at /3 = 45° or -45°, where /3 is the angle between 
Fgtat and the alignment laser field, Faiign [see Fig. Il](a)]. 
The amplitude of the static field is 345 V/cm. We inter- 
pret the S^ ions detected at small radii, near the center 
of the images, as originating from OCS molecules, singly 
ionized by the probe pulse and dissociating into CO and 
S+. By contrast, the S+ ions in the pair of radially and 
angularly localized regions (at the outermost part of the 
images) is interpreted as originating from OCS molecules, 
doubly ionized by the probe pulse and subsequently frag- 
menting into a CO+-S+ ion pair. The recoil of S+ ions 
from this Coulomb explosion channel refiects the direc- 
tion of OCS at the moment of ionization and is thus a 
useful experimental observable to determine the molecu- 
lar alignment and orientation. 

The strong angular confinement of the S+ from 
the Coulomb explosion channel shows that the OCS 
molecules are sharply 1-dimensionally aligned along the 
polarization of the linearly polarized YAG pulse. In addi- 
tion, a pronounced asymmetry of the S+ ions emitted ei- 
ther parallel or anti-parallel to Fgtat, with an excess of S^ 
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FIG. 1: (Color online) (a) Schematic of the velocity map 
imaging spectrometer used to detect ions or electrons. The 
alignment of the molecules, illustrated by the OCS model, is 
determined by the alignment laser polarization, here shown 
for /3 = 45°. The static electric field of the spectrometer, 
pointing from the repeller to the extractor electrode for ion 
detection, breaks the head- for-t ail symmetry by preferentially 
placing the 0-end towards the repeller. When detecting 
electrons the polarity of the electrodes is inverted forcing the 
S-end of the molecules towards the repeller electrode, (b), 
(c) Images of S^ ions for (3 — 45° and -45°, respectively and 
Fstat = 345 V/cm. In (b) 73 % of all S^ ions appear in the 
upper half of the detector. In (c) 28 % of all S^ ions appear in 
the upper half of the detector, (d), (e) Images of S^ ions for 
/3 = 45° and -45°, respectively with the static field increased 
to Fstat = 594 V/cm. In (d) and (e) respectively 80 % 
and 19 % of all S^ ions appear in the upper part of the de- 
tector. The intensity of the probe laser is ~ 5.4 x lO^'^ W/cm^. 



in the upper (lower) region for /3 = 45° (—45°), shows 
that the molecules are oriented with the S-end prefer- 
entially pointing toward the extractor electrode where 
the electrical potential is lowest. These findings are fully 
consistent with recent alignment and orientation studies 
on iodobenzene [29l|30] and 2,6-difluoroiodobenzene [36] 
as well as with former mixed-field orientation studies on 
OCS [37 . In Figs, [ijd) and (e) the static electric field 
is increased to Fgtat = 594 V/cm resulting in a clear 
improvement of the degree of orientation. 

The theoretical treatment of the experimental PAD 
measurements, presented below, requires knowledge of 
the degree of orientation, i.e., the fraction of molecules 



with the S-end pointing towards the repeller. The ion 
imaging measurements occur for beta = 45 degrees rather 
than the degree geometry used in the PAD measure- 
ments. To, nevertheless, provide an estimate of the de- 
gree of orientation we note that Egtat = 345 V/cm in 
the PAD measurements [See Sec. IIIB . This value falls 



in between the value of the effective static field, i.e., 
Egtat along the OCS bond axis, of Figs. 1(b) and 1(c) 
(- cos(45°) X 345 V/cm = 244 V/cm) and Figs. 1(d) 
and 1(e) (- cos(45°) x 595 V/cm = 420 V/cm). In 
the former (latter) case the orientation corresponds to a 
73% (80%) up-to-total ratio. Therefore, the orientation 
in the PAD experiment geometry should be at least 77- 
78 %. The vertical probe geometry applied in Fig. [l] 
does, however, underestimate the degree of orientation, 
because the probe pulse preferentially ionizes (probes) 
the molecules aligned along its vertical polarization axis 
where the static field goes to zero and the molecules are, 
therefore, only weakly oriented. As a consequence, we 
estimate that the orientation in the PAD geometry cor- 
responds to 80 % of the molecules having their 0-end 
toward the detector, see Fig. |2j 



B. PADs from single ionization of OCS 



For the PAD experiments the same experimental setup. 



described in Sec III A is used, but some essential parame- 
ters are changed. The polarization state of the 30 fs probe 
pulses, denoted as F(t) is changed from linear to circu- 
lar and the intensity is lowered to :^ 2.4 x lO"*^^ W/cm^ 
corresponding to a regime where the OCS molecules only 
undergo single ionization with essentially no fragmenta- 
tion. The intensity puts the dynamics in the tunneling 
regime [38 and the circular polarization ensures that no 
recollision of the freed electron with its parent ion occurs. 
Both conditions are important for the interpretation and 
modeling of the observed PADs. Also, the polarization of 
the alignment pulse is changed such that its major axis 
is parallel to the static field axis. Furthermore, to ex- 
tract electrons instead of ions in the PAD measurements 
the polarity of the velocity map imaging spectrometer is 
inverted. Hereby, the OCS molecules are confined along 
the static field axis with the 0-end facing the detector 
(See Fig. [2|. 

The electron images are shown in Fig. [3J With only 
the probe pulse [Figs. |3ja),|3|b)] the electrons emerge in 
a stripe parallel to the (Y ^Z) polarization plane of the 
::^ 2.4 X 10^^ W/cm^ probe pulse. The images are es- 
sentially up-down symmetric and the marginal difference 
between the images obtained with left and right circu- 
larly polarized (LCP and RCP) pulses is due to exper- 
imental imperfections in the purity of the polarization 
state and a weak orientation of the molecules caused by 
the static field alone [364. When the molecules are one- 
dimensionally (ID) aligned along the F-direction, i.e., 
the molecular axis is confined along the F-axis but with 
no preferred direction of the dipole moment, not much 
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FIG. 2: (Color online) Schematic of the experimental setup 
showing an OCS molecule oriented with its permanent dipole 
moment (bold red arrow) pointing in the direction of the static 
electric field. The left circularly-polarized (LCP) probe pulse 
ionizes the molecule and imparts an upward momentum to 
the freed electron resulting in recording on the upper part of 
the detector (see text for details). 



happens and no up-down asymmetry is observed [Figs. 
IsFc) and |3|d)]. When the YAG pulse polarization is 
turned parallel to Fgtat, and the molecules thus become 
ID aligned and oriented, a strong up-down asymmetry 
is observed [Figs, jsje), (sjf)]. The asymmetry reverses 
as the helicity of the probe pulses is flipped. For LCP 
(RCP) probe pulses the number of electrons detected in 
the upper part compared to the total number in the im- 
age is -64% (39%). 

To investigate if the YAG pulse not only induces molec- 
ular alignment and orientation but also influences the 
photoelectron trajectories PADs were measured under 
experimental conditions identical to those used in Figs. 
[3] (e) and [3] (f ) but with the YAG intensity reduced by 
a factor of three. Independent measurements, using S^ 
ion imaging (not shown here) showed that the degrees 
of alignment and orientation remain almost unchanged. 
The resulting electron images are shown in Figs. Is] (g) 
and [3] (h) . They are very similar to those obtained with 
the images obtained at three times higher YAG pulse 
intensity. In particular the up-to-total number of elec- 
trons is r.^66% (38%) in Figs. [3] (g) and [3] (hV, respec- 
tively, which is almost the same as for Figs. [3] (e) and 
[3] (f), strongly indicating that the YAG pulse does not 
cause any significant distortion of the electron trajecto- 
ries. This is corroborated by measurement on benzoni- 
trile where experiments were conducted on molecules at 
higher rotational temperatures and without state selec- 
tion [33^. In that case no up-down asymmetry of the 
photoelectrons is observed even at the highest YAG pulse 
intensity. We conclude that the YAG pulse together with 
the static electric field serve to control the alignment and 
orientation of the molecules but does not otherwise vis- 
ibly influence the ionization process by the probe or the 
subsequent trajectories of the released electrons. 

To investigate the role of the intensity of the probe 
laser pulse, measurements at three different intensities 
were performed. Figure [4] shows the electron images ob- 
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FIG. 3: (Color online) Two-dimensional momentum image of 
electrons produced when a randomly oriented sample of OCS 
molecules are ionized by (a) the LCP probe pulse. The po- 
larization plane of the probe pulse is in the (Y,Z) plane, i.e. 
perpendicular to the image (detector), (b) Same as (a) but 
for a RCP probe pulse, (c) and (d) as (a) and (b) but with 
the OCS molecules aligned along the F-direction by the align- 
ment pulse polarized parallel to the image plane, (e) and (f) 
as (a) and (b) but with the OCS molecules aligned along the 
Z-direction by the alignment pulse polarized perpendicular to 
the image plane, (g) and (h) as (e) and (f) but with the inten- 
sity of the alignment pulse lowered from laiign, yag = 8.4 x 
10^^ W/cm^ to laiign, YAG = 2.8 X 10^^ W/cm^ The inten- 
sity of the 800 nm 30 fs probe pulse is kept at 2.44 x 10^^ 
W/cm^ in all pictures. 



tained with a left circularly polarized pulse. In Figs. 
pFa), [4Fd) and Kg) only the probe pulse is included for 
increasing intensity corresponding to Iprobe = 1.76 x 
IQi"^ W/cm^ in (a), 2.44 x lO^^ W/cm^ in (d) and 
2.83 X IQi^ W/cm^ in (g). In Figs. |4{b),[4];e) andgh) 
the molecules are ID aligned and oriented and again clear 
up/down asymmetries are observed in the photoelectron 
distributions. In Figs. |4]^c),[4]^f) and|4]^i) the correspond- 
ing radial distributions, obtained by angularly integrat- 
ing the images, are given for the upper and lower half of 
Figs. H[b), H[e) andpFh), respectively. It is seen that 
as the probe intensity is increased the number of elec- 
trons detected increases, and they acquire more momen- 
tum, i.e., extend towards the edge of the detector. The 
maximum momentum observed is limited by the size of 
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FIG. 4: (Color online) Two-dimensional momentum 
image of electrons produced when a randomly ori- 
ented sample of OCS molecules is ionized by the 
LCP for increasing intensity of the ionization pulse, 
corresponding to Iprobe — 1.76 x 10^^ W/cm^ 

in (a), Iprobe = 2.44 x 10^^ W/cm^ in (d) and 
Iprobe = 2.83 X 10^^ W/cm' in (g). (b), (e) and (h), 
same as (a), (d) and (g), respectively, but with the OCS 
molecules aligned along the Z-direction by the align- 
ment pulse polarized perpendicular to the image plane. 
In (c), (f) and (i) the radial distributions of the upper 
and lower parts of images (b), (e) and (h) are shown. 

lalign, YAG = 8.4 X 10^^ W/cm" 



the detector. The number of electrons appearing in the 
upper half of the images compared to the lower half is 
almost unchanged for increasing intensity with ~63% in 
Fig. ^h) and -64% in Figs. |^e) and[4];h). These ex- 
perimental observations are compared to the calculated 
results in Sec. IVl 



IV. THEORY 

For comparison with theory, we focus on the intensities 
2.44 X 10^"^ and 2.83 x 10^"^ W/cm^ for the 800 nm, 30 fs 
pulses. For OCS, these laser parameters result in Keldysh 
parameters [38 7 = ujy^/F of 7 = 0.87 and 7 = 0.82 
that are both lower than unity, so it is justified to use 
tunneling theory to describe the photoelectron emission 
process. The existing tunneling models, however, need to 
be modified to correctly describe ionization from a polar 
molecule with large dipole moments and polarizabilities 



such as OCS. 

In a circularly polarized field the electron is driven 
away from the (unrelaxed) cation, which is in contrast 
to the case of a linearly polarized field where rescatter- 
ing and post-ionization interaction are important. This 
fact simplifies the propagation after the initial ionization 
step in circularly polarized fields: it proves sufficient to 
propagate classical equations of motion for the electron 
in the external field ignoring the effect of the molecular 
potential. The full- width at half maximum of the 800 nm 
laser pulse used in the experiment is 30 fs, and accord- 
ingly there are more than 10 cycles within the envelope. 
Therefore it is sufficient to model the laser pulse of the 
experiment by a periodic field with constant amplitude, 
and assume that it is switched off adiabatically in the 
long time limit. We focus on the case of a left circularly 
polarized laser pulse (LCP), and define the electric field 
F(t) of the probe pulse as 



F{t) = Fq sm{ujt)ey + Fq cos{ujt)ez. 



(1) 



In the above equation Fq is the field amplitude, uj is 
the angular frequency and the product 6 = ut is the 
angle between the electric field vector and the positive 
Z-axis. The vector potential corresponding to (IT]) is 
(F(i) = -dtA{t)) 

A(t) = — cos(ujt)ey sm(ujt)ez. (2) 

UJ ^ UJ 

Assuming that there is no influence of the molecular 
potential on the final momentum of the escaping contin- 
uum electron, the emission at time to (angle = ujIq) 
creates an electron with final momentum in the (Y ^Z) 
plane equal to 
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(3) 



From the above considerations it is clear that any 
asymmetry in the electron emission predicted by tunnel- 
ing theory for an oriented polar molecule directly trans- 
lates into the up-down asymmetry (positive-negative Y- 
components of the final electron momentum) observed 
experimentally. Moreover, as it is seen from Eq. (|3|, 
the favoured detection of electrons in the upper half of 
the detector for a LCP field translates into a favoured 
emission from the 0-end of the molecule when the field 
points in negative Z-direction. An obvious candidate re- 
sponsible for the asymmetry in the electron emission is 
the asymmetry of the HOMO in the asymptotic regions 
of large spatial distances (see Fig. p| that enters into 
the tunneling model. To address this question, we first 
briefly review the existing tunneling models. 

The tunneling rate in a static field is governed by the 
exponential exp(-2/^3/(3F)) [39^40^, where F = |F| and 
hz = Y^2/p(0), where lp{0) is the field- free ionization po- 
tential. In the case of an atom in its ground state, the pre- 
exponential factor in the tunneling expression accounts 



for the symmetry of the initial state [40l [41] . An exten- 
sion of the tunnehng theory, fuhy in hne with the atomic 
case, was carried out for the molecules [42]. This molec- 
ular tunneling ionization theory takes the orientation of 
the field relative to the molecular axis into account. The 
rate in molecular tunneling theory [42] for static fields is 



w{F) = 
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(4) 



where B{m') accounts for the orientation of the field with 
respect to the molecular axis. 



B{m') = Y,Ci^Dl,^^{f)Q{l,m). 



(5) 
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FIG. 5: (Color online) Visualization of the degenerate HOMO 
orbitals of OCS, at an iso-density contour value of 0.1. The S- 
end of the molecule is on the left in these images. The orbital 
to the left lies in the polarization plane of the laser field, and 
contributes more to the total ionization yield than the orbital 
to the right, which has a nodal structure in the polarization 
plane and is obtained by 90 degrees rotation of the orbital to 
the left around the molecular axis. The circles with arrows 
pointing counter clockwise illustrate the LCP field. 



In the above equation, D^^, m(-^) ^^ ^^^ Wigner rotation 
matrix element (see, e.g., [43^), for passive rotation of the 
coordinate system through angles F from the molecule- 
fixed frame to the laboratory-fixed frame with the Z-axis 
determined by the direction of the external field. The 
coefficient Q{l,m), given by 



Q{l,m) = {-1) 



i\m\-m)/2 2l + l{l + \m\y. 

V 2 (l-H)!' 



(6) 



is related to the dominant behavior of spherical harmon- 
ics along the field direction [44 , 
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(7) 



In the above equation, 6 and (j) are angular coordinates 
in the spherical coordinate system where the Z-axis is 
directed along the field. Finally, the Ci^rn coefficients are 
related to the asymptotic behaviour of the wavefunction 
of the HOMO, that is 



^(r) 



"^ exp(-A>:r) ^ CimYim{r). 



(8) 
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In Eq. ([8|, the radial part solves the radial Schrodinger 
equation for the electron in the Coulomb field to first 
order in 1/r. The Cim coefficients for both degenerate 
HOMO orbitals of OCS were calculated by projecting 
the orbitals obtained using standard quantum chemistry 
calculation [45 onto the asymptotic form ([8|, yielding 
coefficients with m = ±1 and with / up to / = 5. Since 
the oriented OCS molecule is free to rotate around its 
molecular axis, a combined response from both orthog- 
onal degenerate HOMO orbitals is required. The two 
orbitals are rotated by 90 degrees from each other along 
the molecular axis. We assume one is in the polarization 
plane of the laser pulse and one is perpendicular to it. To 
calculate the combined response from these two orbitals 



the associated angle-dependent tunneling rates of both 
orbitals are added incoherently since the molecule is in a 
mixed state with respect to the degenerate orbitals. 

Direct application of the molecular tunneling theory 
sketched above, however, gives an opposite emission de- 
pendence to the one required to describe the asymmetry 
observed in the experiment. Namely, the B coefficients 
from ([5| are such that the emission from the S-end is 
favoured with respect to the emission from the 0-end. 
This reflects the simple fact that the HOMO orbitals of 
the OCS are such that the wavefunction is predominantly 
located toward the S-end of the molecule (see Fig. [s]). As 
we will show below, however, to describe the tunneling 
ionization process from a polar molecule it is essential 
to take into account the angle-dependent shifts of the 
ionization potential, induced by the polar system. Un- 
like the 5-coeflicients that depend on the geometry of 
HOMO and influence the pre-exponential factor in the 
tunneling rate of Eq. Q, the Stark shifts affect both the 
pre-exponential factor and the argument of the exponen- 
tial in the tunneling rate and decide from which end of 
the molecule the preferred emission occurs. 

Our modiflcation of the tunneling theory is building 
on the fact that in electric flelds, due to the molecule's 
polarity and its polarizability, hyperpolarizability, etc., 
the energy levels of the molecule shift. These shifts are 
negligible for small nonpolar molecules and atoms. In the 
case of molecules with large dipole moments and polariz- 
abilities the influence of these Stark shifts cannot be ne- 
glected and in the present quasistatic limit they must be 
included. The unrelaxed cation of the molecule in ques- 
tion is even more tightly bound than the neutral molecule 
and the characteristic timescale for the electronic motion 
is shorter. Hence, if for the molecule the fleld can be 
regarded as static, then for the cation the fleld is static 
as well. The total energy of a molecule (M) and its unre- 
laxed cation (I) E^^^ {¥) in a static fleld F, up to second 



order in field strength is given by (see, e.g. |1S]) 

^M//(p) ^ ^M//(Q) _ ^M/i . p _ IpTa^/^F, (9) 

where ii^^^ is the dipole moment, oJ^I^ is the polariz- 
abihty tensor, £^^/^(0) is the field- free total energy of 
the system. The next term in the expansion of the total 
energy as a function of field strength involves the hyper- 
polarizability. In OCS, the contribution of the hyper- 
polarizability is negligible at the intensities used in the 
present experiment, so the Stark shift is due to the dipole 
moment and the polarizability only. The Stark shift due 
to the polarizability is larger than the Stark shift due to 
the permanent dipole moment. What is important, how- 
ever, is the difference between the total energy of the 
molecule and the positive ion, i.e., the ionization poten- 
tial. Since the molecule and the ion do not have identical 
permanent dipole moments and polarizabilities, the ion- 
ization potential Ip = E^ — E^ becomes 

/^(F) = /^(O) + A/i • F + ^F^AaF, (10) 

where F^ is the transpose of the field vector and 

Afi = fi^-fi' Aa = a^-a'. (11) 

In the above formulation, we assumed that the electron 
in highest occupied molecular orbital was promoted to 
the continuum. In this case the change of the ionization 
potential as a function of field strength can be referred 
as to as the Stark shift of the HOMO orbital, with the 
corresponding dipole moment and polarizability. Such 
modifications of the ionization potential of orbitals other 
than HOMO can be calculated as well. Note that here we 
do not make any distinction between the dipole moments 
and polarizabilities of (possibly) different ionic products 
obtained through different ionization channels. We sim- 
ply refer to the properties of the cation in the unrelaxed 
geometry of the neutral molecule. 

Equations (10) and (11) show very explicitly that the 



ionization potential depends not only on the magnitude 
of the electric field vector F but also on the angles of 
the field orientation with respect to the principal polar- 
izability axes and the permanent dipole moment of the 
molecule. This is the essential ingredient that enters into 
our modification of the tunneling theory. Note that static 
Stark shifts were considered earlier in discussion of dis- 
sociation [13 [48] 

The modification of the ionization potential due 
to Stark shifts results in the tunneling exponential 
exp(— 2ai:^(F)/3F), where the factor 



k(F) = J2Ip{F) 



(12) 



initial orbital is not affected by the polarization of the 
molecule. This is in general not true, especially for field 
strengths at which the polarizability term gives much 
larger contribution to the Stark-shifted ionization poten- 
tial of Eq. ( 10 ) than the permanent dipole moment term. 



dependens on the angles of field orientation with respect 
to the molecular axis. One can now in principle take the 
tunneling rate of Eq. (pi) and replace k, with k,{F) every- 
where. However by doing so, one has assumed that the 



In these cases, the polarizability can modify the initial 
molecular orbital so that the Cim coefficients in Eq. (Isl) 
become a function of the field strength. It is in general 
very hard to account analytically for such strong modifi- 
cations of the initial orbital. However, as we show below 
by direct comparison to the experiment, in the case of 
OCS and in the limit of large fields it suffices to sim- 
plify the situation by disregarding the modifications of 
the HOMO and assuming that the angle-dependence of 
the tunnel emission occurs only due to the action of the 
Stark shift. Namely, we assume that the inner structure 
of the orbital is modified so much by the polarization re- 
sponse so that its asymptotic properties would be equiv- 
alent with respect to the axis defined by F for each orien- 
tation of the field with respect to the molecular axis. In 
other words, if the Stark shift would not be present, there 
would be no orientation dependent emission, i.e., the tun- 
neling probability would be equal for all field orientations 
with respect to the molecular axis. Having this in mind, 
we can model this behaviour by taking an atomic s-like 
state as an initial state in the tunneling model. This 
radical model will generally be better for large intensi- 
ties, for systems with large polarizabilities and in cases 
where the initial orbital is not such that the polarization 
plane of the circularly polarized laser field lies entirely in 
the nodal plane of the orbital [33] . 

The intensity of the laser pulse used in the experiment 
is very large so that for atoms with the same binding 
energy the ionization would occur over the barrier. To 
examine whether that holds for the case of the OCS 
molecule one requires a single-active-electron potential 
corresponding to the HOMO of OCS. Following the ap- 
proach given in Refs. [19l[49], applied sucessfully to CO2 
and other linear molecules, we have built a single- active- 
electron potential for OCS. 

Using the single-active-electron potential for OCS, we 
have verified that at the peak intensities of the experi- 
ment, the ionization occurs over the barrier at all angles 
of orientation of the field with respect to the molecu- 
lar axis and both with or without inclusion of the Stark 
shift. The saddle points of the potential at the experi- 
mental peak intensity of 2.44x10^^ W/cm^ are shown in 
Fig. [6| It can be seen that the saddle points occur rela- 
tively far away from the center-of-mass coordinate so the 
infiuence of the potential to the outgoing electron, born 
at the saddle points is very small. Hence the asymme- 
try of the molecular potential does not play a role in the 
present case. In addition, due to the large polarizabil- 
ity of the parent ion (see Table [H] in Appendix A), the 
induced dipole of the cation is very large and the orienta- 
tion is such that it shields the electron from the infiuence 
of the attractive Coulomb potential, thereby decreasing 
the effect of the long-range Coulomb potential. 
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FIG. 6: (Color online) The saddle points of the OCS potential 
at the intensity of 2.44 x 10^^ W/cm^ in the polarization (Y-Z) 
plane. The dots are the saddle points of the effective potential 
made up of the molecular single-active electron potential and 
the external field. The values of the saddle points are taken at 
5 degree steps with respect to the z-axis, and are connected 
with the dotted lines to guide the eyes. The origin is the 
center of mass of the molecule. The S-end of the molecule is 
to the left of this figure, in accordance with the geometry of 
Fig. [2] The explicit nuclear positions are as follows. S (-2 
a.u.), C (0.99 a.u.), and O (3.2 a.u.) (see Appendix A). 



At over-the-barrier intensities, the tunneling rate given 
by the tunneling theory overestimates the probability 
of ionization [50l |5j. The saturation of the tunneling 
rate can be included by calculation of the exact, com- 
plex eigenenergies in a static field [511 [52], however this 
approach is limited only to very simple systems. In this 
study, we will adopt an ad-hoc and simple exponential 
factor, given in Ref. [50 , to account for the over-the- 
barrier saturation of the tunneling rate. The saturation 
factor due to over-the-barrier emission reads. 



V^(F) =^(F)exp -6 



i^^{F)J V^^(F) 



(13) 



where we have additionally included the Stark shifts in 
K according to Eq. ( [T2| ) since it occurs as an argument 
in the exponential. In the above equation, w{F) is the 
tunneling rate from Eq. Q, where Ci^rn — ^^o^m,o cor- 
responds to an s-state and k has been replaced by /^(F). 
In summary, in the tunneling model we use to describe 
the momentum distributions for the OCS molecule, the 
tunneling rate is calculated as 



W(F) 



1 



2a^(F) 
2/^(F)^"^ V ^ 
2 



3\ ^(F)" 



exp 



2k{FY 



3F 



X exp —6 



M^f W{F)' 



(14) 



COMPARISON WITH THE EXPERIMENT 
OF STRONG-FIELD IONIZATION OF 
ID-ALIGNED OCS 



We consider the experiment on strong-field ionization 
of ID oriented OCS by a circularly polarized field in the 
(y, Z) plane. The molecular axis is along the Z-axis, and 
it is perpendicular to a detector that lies in the (X, Y) 
plane (see Fig. [2|. 

We assume the electric field of a LCP pulse as given 
by Eq. ([l]) and include the static Stark shift of the active 
HOMO through the shifts of the molecule and the unre- 
laxed cation as given in Eqs. (lO)-(ll). This approach 



leads to a modification of the ionization potential as a 
function of the angle of the direction of the field F with 
respect to the Z-axis (denoted as ^ = F = ujt) as 



lp{0) 



lp{0) 

1 

+ 2 

-*■ -72/ M 



(/-//^)Focos(0) 



F^ {{a^ 



zz 



-M^ 



XX 



^zz) 



^xx)' 



yo^xx 



(15) 



with hi(Y) given by Eq. ([T2|). 



Note that for the angle of the field, the emission oc- 
curs in the direction 6> + 7r. The above equation gives the 
modification of the ionization potential as a function of 
the angle between the instantaneous direction of the field 
and the permanent dipole of the molecule, and together 
with Eq. (14) it is the main theoretical input in the in- 
terpretation of the experiment. The modification of the 
ionization potential is such that when the field vector and 
the permanent dipole moment of the molecule are paral- 
lel, the ionization potential is minimal so the tunneling 
probability reaches maximum. In that case the emis- 
sion occurs opposite to the field direction, that is, from 
the O-end of the molecule. On the other hand, when 
the instantaneous field points from the S- to the O-end 
of the molecule, i.e., the field and the permanent dipole 
moment are antiparallel, the ionization potential is max- 
imal and, conversely, the ionization probability minimal. 
The emission then occurs from the S-end of the molecule. 
In the simplest model, the final momenta in the {Y- 
Z) plane are given by Eq. p|. Since there is larger 
probability of tunneling from the O-end of the molecule, 
for the LCP, there would be a larger probability that 
electrons with py > appear, and for RCP pulse the 
situation would be opposite, which is in accord with the 
experiment. Scanning through instants of time within 
one cycle of the field, all possible final momenta py and 
pz in the plane px = are reached and the associated 
ionization probability calculated from Eq. (14), where 
|F| = Fq and the orientation of F is given by defined 
in Eq. ([3]). On the other hand, the transverse, px com- 
ponent of the momentum cannot be changed by the ex- 
ternal electromagnetic field and it is obtained from the 
well-known expression [53] for the momentum distribu- 
tion of the transverse momenta of the tunneled electron 
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FIG. 7: (Color online) Momentum ditributions obtained by 
the present model for peak intensities of (a) 2.44 x 10^^ and 
(b) 2.83 X 10^^ W/cm^ of the 800 nm, 30 fs probe laser pulse. 
Compare with the experiments in Figs. H] (e) and (h). 



at its birth, i.e., 



W{px) ~ exp I -^^ p^ 



(16) 



Note that in actual calculations for the transverse distri- 
bution the approximation lp{0) ^ ^p(O) has been used 
since the inclusion of Stark shifts in this degree of free- 
dom does not alter the results significantly and has no 
influence on the observed asymmetry in the experiment. 
In the experiment, the measured quantity on the de- 
tector is 



W{px,Py) = 
where VK(p) 



dpzW{p), 






dpxdpydpz 



(17) 



and Pion is the total ionization probability. In the sim- 
pleman model, by inserting Eq. (15) into (14), the 



tunneling rate W{6) can be obtained, and from there, 
W{^^PYiPz) = W{9)^ where py and pz dependency on 
6 is taken from ([3|. This momentum distribution is then 
integrated over pz to obtain k;(0,py). Having in mind 
Eq. (16), the momentum distribution (17) is obtained as 



W{px,Py) = W(0,py)exp 



72^(0) o 



(18) 



As discussed in section III A[ the orientation of the 
OCS molecules in the experiment is not perfect. In fact, 
80% of the molecules are oriented in the desired orienta- 
tion (0-end towards detector) and 20% oppositely. These 
20% of the molecules actively participate in the formation 
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FIG. 8: (Color online) Comparison of experimental (full lines) 
and theoretical (dashed lines) radial distributions for peak 
intensities of (a) 2.44 x 10^^ and (b) 2.83 x 10^^ W/cm^ of 
the 800nm, 30 fs probe laser pulse. The upper curves are 
the radial distributions in the upper (jpy > 0) and the lower 
curves are radial distribution in the py < half plane, see 
Eqs. ([22l-([23l. 



of the experimentally obtained momentum distribution 
hence this is an effect which must be taken into account 
in the theoretical model. The corresponding momentum 
distribution with this effect taken into account is readily 
obtained from (18) as 



WpniPx.Py) = 0.8W{px,Py) + 0.2W{px, -pv)- (19) 

The final effect which must be taken into account to 
reproduce the experiment is the volume effect for tightly- 
focused laser beams. The consideration of this effect 
requires the calculation of momentum distributions of 
type (19) for different intensities and then using a vol- 



ume function to weight the momentum distributions at a 
particular intensity [54J . After inclusion of the volume ef- 
fect, we obtain the momentum distributions given in Fig. 
nfor the experimental peak intensities of 2.44 x 10^^ VV^)] 
and 2.83 x 10^^ W/cm^ Rt))]. These momentum distri- 
butions visually resemble very much the corresponding 
experimental momentum distributions of Figs. [4] (e) and 
(h). As with their experimental counterparts, the calcu- 
lated momentum distributions in Figs, [t] (a) and (b) are 
very similar, reflecting that the momentum distributions 
are only weakly dependent on the laser intensity. 

A more detailed comparison between the experimental 
[Figs [4] (e) and (h)] and theoretical momentum distribu- 
tions [Figs, [t] (a) and (b)], reveals larger F-components 
of the final momenta in the theoretical momentum distri- 
butions, see also Fig[8J This is due to two reasons, which 
are not included in the theory presented here. First, 
the target OCS molecules are not ideally aligned with 
(cos^ 0) ~ 0.9. Hence, there is a non-vansihing prob- 
ability that the angle between the molecular axis and 
the space-fixed Z-axis is nonzero. If we recall the post- 
ionization dynamics discussed in the previous section, the 
continuum electrons with the largest F-component of the 
final momentum escape into the continuum at times when 
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FIG. 9: (Color online) Comparison of experimental (full lines) 
and theoretical (dashed lines) angular distributions [Eq. (|21|) 



for peak intensities of (a) 2.44 x 10^^ and (b) 2.83 x 10^^ 
W/cm^ of the 800 nm, 30 fs probe laser pulse. 



FIG. 10: (Color online) Comparison of experimental (dots) 
and theoretical (dashed lines) asymmetry A'^ [Eq. ( |24[ )] as a 
function of radial momentum for peak intensity of (a) 2.44 x 
10^^ and (b) 2.83 x 10^^ W/cm^ of the 800 nm, 30 fs probe 
laser pulse. 



the electric field vector F(t) is parallel or antiparallel with 
the Z-axis [see Eq. ([3|]. On the other hand, emission 
from the 0-end of the molecule is preferred and due to 
the nonzero value of the angle between the space-fixed 
Z-axis and the molecular axis, the absolute value of the 
F-components of the final momenta would be smaller 
than the maximal value, i.e. \py\ < Fq/uj, so the Y- 
components of the momenta are lowered. The second 
reason for small \py \ in the experimental momentum dis- 
tributions is the interaction of the outgoing electron with 
its parent ion, which is beyond the simpleman model. 
Although in the case of a circularly-polarized laser field 
such post-ionization interaction is limited, the long-range 
Coulomb part of the potential will act to decrease the 
magnitude of the final momenta due to the attractive 
forces at times immediately after ionization. 

The above effects, however, have no influence on the 
most prominent feature observed in this experiment - 
the up-down asymmetry in the momentum distributions. 
This asymmetry perpendicular to the permanent dipole 
moment is entirely described by including the static Stark 
shifts in the tunneling model. The up/total asymmetry 
A+ is defined as 



A+ 



)y>o 



PZ^VP.o.. where 



ton I ^ ton 

^oo /»oo 



p — 



/CO />00 

dpx / dpvWpnipx.PY) and (20) 
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/oo />oo 

dpx / dpY Wpn {px ,Py)' 
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We have calculated A^ = 0.651 (experimental value 0.64) 
at the intensity of 2.44 x 10^^ W/cm^ and A+ = 0.649 
(experimental value 0.64) at the intensity of 2.83 x 10^^ 
W/cm^ . We conclude that the theoretical and the exper- 
imental asymmetry agree well. 

We turn to the comparison of the differential quantities 
that can be derived from the momentum ditributions. 



Transforming the momentum distribution of Eq. ( 19 ) into 
polar coordinates px = Pp cos((/)) and py = Pp sin(^, one 



obtains the angular distribution 



/>oo 

Jo 



(21) 



shown in Fig. [9] and the radial distribution in the upper 
half {pY > 0) 



WP->'{p, 



-PpI d(l)W{pp,^) 





and in the lower half of the plane {py < 0) 



WP^<\p,)=pJ d^W{p, 



(22) 



(23) 



shown in Fig. [8) 

The results for the angular distribution in Fig. [9] are in 
excellent agreement with the experiment. In particular 
the peak ratio of the peaks in the angular distribution 
for 90 and 270 degrees is reproduced by the theory. The 
width of the two peaks in the experimental radial dis- 
tributions is slightly larger than the theoretical peaks. 
This is due to the non-perfect alignment of the molecule 
along the Z-axis discussed above. Moreover, one has to 
consider Coulomb focusing: the electron wavepacket is 
attracted to the (F, Z) plane, which results in a narrower 
width of the angular distributions. Thus the nonperfect 
alignment and the effect of the Coulomb focusing have an 
opposite effect on the angular distributions, almost can- 
celling each other and resulting in an overall very good 
agreement between theory and experiment. On the other 
hand, as discussed above, the non-perfect alignment and 
the long-range interaction on the outgoing electron both 
have a tendency to decrease the final momenta. This 
is evident from the comparison of the experimental and 
theoretical radial distributions presented in Fig. [9] The 
theoretical model overestimates the F-components of fi- 
nal momenta. We note that the sudden cut-off of the 
experimental radial distributions is due to the maximum 
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momentum that can be recorded on the detector, and 
not because of some physical effect. 

The final quantity that has been calculated from the 
experimental data is the up/total asymmetry as a func- 
tion of the radial momentum A+(pp), defined as 



Jo Jo 



,0). (24) 



Experimental and theoretical results for this quantity is 



presented in Fig. 10 For both intensities, the results for 
the differential asymmetry A'^{pp) agree very well in the 
region of smaller radial momenta. Namely, the differen- 
tial asymmetry rises gradually from the value of 0.5 (no 
asymmetry) to around 0.7 at the peak. Very small radial 
momenta correspond to ionization when the electric field 
vector points in the direction of the positive or negative 
y-axis (electric field vector perpendicular to the dipole 
moment) with equal probability of ionization in both di- 
rections [put 6 = 7r/2 in Eq. (15)], therefore no asymme- 
try. As the radial momentum pp increases, \py\ increases 
and I cos(^)| of Eq. (fTsl) also increases, resulting in larger 
asymmetry. The theoretical model captures these fea- 
tures in the experimental data. The small discrepancy of 
experimental and theoretical results for A'^{pp) at small 
and intermediate Pp values is again due to the overesti- 
mation of the Y components of the final momenta by the 
theory, discussed above. At the largest Pp values the de- 
crease of the asymmetry in the experimental curves is not 
reproduced by theory. This could be due to nonperfect 
alignment or to focal volume effects, i.e., the asymme- 
try at 0.6 comes not only from the contribution of the 
tunneling slightly off the O end of the molecule at the 
peak intensity but also from tunneling exactly from the 
O end for lower peak intensities contributing to the total 
signal. Finally, the discrepancy could arise because the 
model neglects any orientation-dependent tunneling that 
has its origin in the initial state. 



VI. CONCLUSIONS 



In DCS and for the present set of laser parameters, the 
electron enters the continuum at such large distances that 
the asymmetry associated with the ionization potential 
in the initial tunneling process is sufficient to explain 
the experimental findings. The asymmetry associated 
with the molecular potential does not play a role, since 
the leading asymmetric dipole term is suppressed at the 
distances in question [see Fig. 6]. 

It is our goal to extend the present techniques to study 
time-resolved electron dynamics, for example, to moni- 
tor changes in the MFPADs during a photochemical re- 
action. In such a process the nuclear motion could be 
sufficiently slow that changes in the MFPADs will oc- 
cur at the femtosecond timescale [55] , and changes could 
then be recorded by firing a short femtosecond pulse. In 
the tunneling regime, asymmetries in the MFPAD will 
be exponentially sensitive to the ionization potential and 
hence to changes in the dipole moments and polarizabil- 
ities. Changes in the MFPADs will consequently link 
directly to the instantaneous values of these quantities. 
There might be changes in the MFPADs for reasons other 
than the ones connected with the dipoles and polariz- 
abilities. We mentioned asymmetries due to the exact 
form of the molecular potential above. Another possibil- 
ity is the temporal formation and changes of nodal sur- 
faces. It is encouraging that signatures of nodal surfaces 
is clearly detectable with the present technique, as shown 
for benzonitrile [3TJ |33] , and as discussed theoretically in 
Ref. [56]. 
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We have studied photoionization of aligned and ori- 
ented OCS molecules. The prepared gas sample was ion- 
ized by intense near-infrared femtosecond laser pulses. 
The detection of ionic fragments was used to investigate 
the degree of alignment and orientation. The sample 
was then used to investigate molecular frame photoelec- 
tron angular distributions (MFPADs) in the tunneling 
and over- 1 he-barrier regimes. Strong asymmetries in the 
distributions were observed and explained in terms of 
a modified tunneling theory. The circularly polarized 
field steers the electron away, and minimizes rescatter- 
ing, but the presence of a permanent dipole moment and 
a large polarizability of the active orbital, the HOMO 
in the present case, means that the effective ionization 
potential shifts depending on the instantaneous magni- 
tude and direction of the external field with respect to 
the molecular axis. 



APPENDIX A: Molecular properties of OCS 

The molecular properties of OCS were obtained for the 
experimental geometry [53 based on the Hartree-Fock 
wavefunction [45 , with the molecule oriented along the 
Z-axis and the 0-end pointing towards the detector (see 
Fig. [2]). In the coordinate system, defined in Fig. [2I 
where the center of mass of OCS is the origin, the Z- 
coordinates of the atomic centers are -2 a.u. for S, 0.99 
a.u. for C, and 3.2 a.u. for the O atom. The first ion- 
ization potential (/p), computed as the positive energy 
of the highest occupied molecular orbital (HOMO), is 
11.4 eV (0.42 a.u.), in agreement with the experimen- 
tal value of 11.2 eV [58]. The OCS molecule has an 
asymmetric charge distribution with a dipole moment of 
0.71 D (1 a.u.=2.54 D) pointing towards the S-end [59^. 
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TABLE I: The dipole moment of the OCS molecule and the 
cation (in the geometry of the neutral molecule and both 
pointing from the 0-end to the S-end) , given in units of Debye 
(1 au=2.54 Debye). 

State TZV" ACCT^ PC" APC^ Sadlej" 



OCS 
OCS^ 



0.52 
1.83 



0.72 
2.69 



0.71 0.88 
2.00 2.20 



0.65 
1.61 



^Valence triple-^" basis set |60| . 

^Dunning triple-^" basis set [6j. 

^Jensen polarization-consistent basis set |62| . 

^Jensen polarization-consistent basis set with diffuse basis func- 
tions [63| . 
^Sadlej pVTZ basis set |64j- 



TABLE IL The non-zero components of the polarizability and 
hyperpolarizability for the OCS molecule and the cation (in 
the geometry of the neutral molecule), given in atomic units. 

State axx o^yy ctzz Pxxz I^yyz P>zzz 
OCS 26.15 26.15 50.72 -45.92 -45.92 -12.85 
OCS+ 19.06 18.73 44.09 -17.55 -20.17 18.49 



The OCS molecule has two degenerate HOMO orbitals 
shown in Fig. [5] By expanding the total wave function 
in a linear combination of atomic orbitals ^i, centered at 
each nucleus ^ = ^^ Qc/)^, we obtain an estimate of the 
electron population at each center as (q)^. This popula- 
tion analysis indicates that 75% of the HOMO electron 
density is localized on the S-end, 15% on the 0-end and 
only 10% on the C atom. The orbital lying just below the 
HOMO in energy has a significantly higher ionization po- 



tential (17.1 eV from Hartree-Fock calculations [IS]), and 
since we are in the tunneling regime with an exponential 
sensitivity to the ionization potential its contribution to 
the ionization dynamics is expected to be negligible. 

In the development of the tunneling theory, we need 
the dipole moments, polarizabilities and hyperpolariz- 
abilities of the OCS molecule. These were first computed 
using the Hartree-Fock wavefunction in conjunction with 
the aug-cc-pVTZ [GT basis set. Comparisons to results 
reported for the neutral molecule in the Computational 
Chemistry Comparison and Benchmark DataBase [ 65] , 
suggest that while the polarizabilities and hyperpolariz- 
abilities are reasonable, the dipole moments should be 
improved. We carried out calculations using the MP2 
level of theory with different basis sets, and the resul- 
tant dipole moments are shown in Table |T| For the neu- 
tral molecule, the computed dipole moments are in fair 
agreement with the experimental value (0.71 D [59 ). For 
the cation, on the other hand, we obtained a wider range 
of values depending on the basis set. A simple classical 
model based on the redistribution of atomic charges due 
to removal of the HOMO electron, without allowing the 
remaining electrons to relax, predicts \iz = -2.2 D for the 
cation. This is in very good agreement with the results 
based on the Jensen's polarization-consistent basis sets 
( -2.0 Debye without diffuse functions [62 vs. -2.2 De- 
bye with diffuse functions [63 ). In the tunneling model 
presented in Section [IV| we use the dipole moment for 
OCS+ obtained via the simple charge model. The molec- 
ular polarizabilities and hyperpolarizabilities are shown 
in Table El 
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